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We calculate the time independent four-point function in high temperature (T) QCD and obtain 
the leading momentum dependent terms. Furthermore, we relate these derivative interactions to 
derivative terms in a recently proposed finite T effective action based on the SU(3) Wilson Line 
and its trace, the Polyakov Loop. By this procedure we thus obtain a perturbative matching at 
finite T between QCD and the effective model. In particular, we calculate the leading perturbative 
" l ] QCD-correction to the kinetic term for the Polyakov Loop. 
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(N . I. INTRODUCTION 

At high temperatures QCD is expected to be found in a new phase, the quark-gluon plasma. While the thermal 
excitations are hadrons and glueballs at low T, the degrees of freedom in the plasma phase are the quarks and gluons. 
This new state of matter is believed to have existed during the first microseconds after the Big Bang, and much of 
the recent interest stems from the fact such conditions may be produced in heavy- ion collisions. Already the results 
from CERN-SPS seem to hint in that direction, and the experiments at higher energies at BNL-RHIC have provided 
^vq i a wealth of new interesting results after its first year of running Q . To understand and interpret the experimental 
| signatures in terms of the evolution of the initial stage after a heavy-ion collision is clearly a very challenging theoretical 
i~ I . task. 

The most convincing theoretical results that a drastic change in the degrees of freedom takes place at a certain T 
come from lattice studies. In the pure glue theory, there is a phase transition between the confined and deconfined 
phases at a critical temperature T c ~ 270 MeV When massless quarks are added, there is similarly a phase 
transition to a chirally symmetric phase at T c ~ 155 — 175 MeV where the precise value depends on the number 
of flavors. Recent lattice simulations also suggest that the chiral transition is simultaneous with the deconfining one 
At physical quark masses the situation is not completely clear, in the sense that there may not be a true phase 
transition but only a rapid cross-over ||. Nevertheless, lattice simulations have shown that the pressure, divided by 
the ideal gas result and plotted against T/T c , is almost independent of the number of flavors Q. 

Due to asymptotic freedom, the quark-gluon plasma behaves as an ideal gas at asymptotically high temperatures. 
Up to corrections of the order of 20 percent, this behavior holds down to temperatures T ~ 3T C , even though each 
higher order term in a straightforward perturbative expansion gives widely different contributions in this temperature 
regime wW- Instead, at T < 5T C one needs a resummed effective theory in terms of quasi particles, the HTL effective 
action B. Such an effective description correctly reproduces thermodynamic quantities like the pressure, as measured 
^ ■ by the lattice, down to approximately 2T C ]io| , [Tl[ |. 

Despite this progress, it is of course highly desirable to actually have an analytical description at T ~ T c , close to 
the critical temperature. Since the QCD coupling constant g ~ 2.5 at T c (using a renormalization scale /i = 2ttT), 
one is presumably forced to consider effective models that go beyond the fundamental QCD Lagrangian. In a recent 
paper p2|, such an effective theory was constructed in terms of the thermal Wilson Line L, 



o 



43 
Oh 



43 



L = V exp 



■13 

ig I drAo (x, r) 



(1) 



where V denotes path ordering, (3 is the inverse temperature and Aq = AqT 11 the time component of the gluon 
field, with T a the generators of the fundamental SU(3) representation, a — 1, ... ,8. The trace of the Wilson Line is 
proportional to the Polyakov Loop /, I — (l/3)TrL. In the pure Yang-Mills theory, the Polyakov Loop is an order 
parameter for a global Z(3) symmetry separating the confined and deconfined phases, with (/) ^ ((I) =0) above 
(below) the phase transition |Dj]. When dynamical quarks are introduced, / ceases to be an order parameter in the 
strict sense, but the susceptibility of I still peaks strongly at T c p|,|l4j|. 

In the effective theory |12| , the pressure of the quark-gluon plasma at T > T c is completely due to the condensate 
of I, and below T c , where(Z) = 0, the pressure vanishes. Moreover, the effective potential V(l) changes extremely 
rapidly around T c . Hence, as the system cools it may find itself trapped at the wrong value of (l). By coupling the 
effective field I to hadronic degrees of freedom, e.g. the pions, hadrons can be produced as I evolves from (I) ^ 
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and subsequently oscillates around (I) = 0. This scenario is somewhat reminiscent of reheating after inflation [fLq 



and much attention has lately been paid to that aspect of the model 1 16 IT] . Remarkably, many qualitative features 
observed at RHIC are in accordance with the model predictions. 

However, when it comes to questions related to the change of the expectation value of I and particle production 
around T c , one has to take into account the variation of / in space-time. In this paper we address the question of 
radiative corrections to the spatial variation, by considering the leading one-loop QCD contribution to the spatial 
derivatives of I. Previous work K,^7) took into account only the classical kinetic term in the (Euclidean) effective 
action, T(l) — (l/2){|<9 t Z| 2 + |<3j£p} + V(l). While such an approach certainly is justified at these preliminary stages, 
it is important to estimate how much the radiative QCD-effects can affect T(l) around T c , where the QCD coupling 
constant becomes large. As for the parameters in the potential V(l), they can be fitted by comparing to QCD lattice 
results and so are well defined at all T. The kinetic term, on the other hand, has to be matched to perturbatively 
calculated terms in QCD, and could therefore receive large radiative corrections. The magnitude of the first one-loop 
QCD correction can then hopefully serve as a guideline to the importance of loop effects, and indicate how reliable 
the above form of T(l) is at T c . We want to stress that the radiative corrections to be discussed come from QCD, and 
not from fluctuations in T(l). 

To study the correction to the kinetic term \dil\ 2 , we first consider the one-loop induced quartic terms in QCD, that 
contain four powers of the external field Aq and two powers of the external momenta. These terms contribute to the 
high T, dimensionally reduced QCD effective action T(^4o) j[8|,[l9), and apart from providing a correction to T{Aq) 
they can also be related to the kinetic term in T(l). As a byproduct we obtain some additional derivative interactions 
in T{1). 

The paper is organized as follows. In the next section, we give the perturbative QCD calculation that corresponds 
to the leading derivative interactions in T(A ). In Sec. Ill we make the actual matching from an effective theory in 
terms of Aq to the one in i, and discuss the validity of the results. We end with our conclusions and an outlook. Our 
conventions and some technical details can be found in the appendix. 

II. PERTURBATIVE CALCULATION OF THE FOUR-POINT FUNCTION 

In the high temperature regime, long distance phenomena (i.e. \x\ ^ (3) are dominated by the static sector of QCD. 
At high T it therefore makes sense to use dimensional reduction and integrate out all the nonstatic modes in the 
theory M]. With only the static modes left, the full QCD Lagrangian is reduced to a three-dimensional theory. In 
principle the integrating-out procedure gives rise to an infinite number of interaction terms, but higher dimensional 
operators become more suppressed by powers of the QCD couplin g co nstant g and/or T. In full QCD, the following 
terms in the resulting effective action T(Aq) have been calculated |l8|]l9f l, 



T(A ) = / d 6 x 



iurfg +Tr [A,A)][A, M + g 2 T 2 (l + Zf\ Tr A% + 9 \j f) A 2 f 



(2) 



where i, j = 1, 2, 3, = F^T a = (d l A a j - d 3 Af - gf ahc A\A<j)T a , D d = di + igA t and A. t = A°;T a . The next term in 
T(Aq) contains two derivatives and four powers of Aq, and corresponds to the following part in the Euclidean effective 
action for Aq(x), 

4Vo) = ^ S(3) ^ + ••• + **) [-irgggfr*;, • • • , £)] A5(&)A>&)AZ{k 3 )A*(ki) , (3) 

where rg§QQ is the four-point function of order 0((3 2 k 2 ), obtained by integrating out all the non-static modes. Such 
higher dimensional terms have been calculated in the pure Yang-Mills theory using the background field method |20| , 
as well as in QED p9| , but not in QCD with quarks. In this paper we will use a diagrammatic approach to the 
four-point function. 

Perturbatively, the four-point function receives contributions from the diagrams shown in Fig. 1, together with the 
additional permutations of the external legs. There are five permutations adding to the graphs (a), (b) and (d), and 
two to the diagrams (c) and (e), where (e) has a symmetry factor 1/2. Even though all diagrams are superficially 
logarithmically divergent, it is well known that both the fermion diagram and the sum of the pure Yang-Mills diagrams 
are ultra-violet finite. We will therefore only give explicit results for the finite T part of these diagrams, where we use 
the imaginary time formalism p2| combined with the particular technique described in the appendix. 
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FIG. 1. Contributions to Foooo-' (a) fermions, (b) ghosts and (c), (d), (e) gluons. 



A. The fermion contribution 



Consider first the contribution from the Nt massless quarks, with the particular ordering of momenta shown in Fig. 
1(a). Using the conservation of momentum, k^ — — {k\ + k 2 + k 3 ), with fej = fcj, we find 



re+ioc j 



— (1 - 2 1 - 2 )r(z)^(z) cos[7rz/2]r z x 

e— loo ^ 7rl 



«i (Q)1 - (^j t3tc(t^pt-) y o ^ y y o y 

/r z/2 r(z/2) ||r((i-z)/2) - er((3-z)/2) + Hr((5-z)/2)| + /r (2+z)/2 r((2+z)/2) 
/ 3 r((i-z)/2) - H/ 2 r(( 3 -z)/2)\ + ^/ 4 /r (4+z)/2 r((i-z)/2)r((4+z)/2) 



(4) 

where £(z) is the Riemann Zeta-function and /i = /i(t/i, 2/2> 2/3) &i> • • • > fe); given explicitly in the appendix, are 
functions of the Feynman parameters yk and the external momenta kj. In the high temperature limit, where \ki\ <C T, 
we can use the residue theorem to evaluate the z-integral by closing the contour on the left side in the complex z-plane. 
Although there is seemingly a logarithmic dependence on T from a double-pole at z — 0, coming from the product 
r(z)r(z/2), the coefficient is actually proportional to [(5/2)r(l/2) - 6r(3/2) + (2/3)T(5/2)] = 0. This is in accordance 
with the fact that the fermion loop does not have any logarithmic enhancements [ |23fl . 
For the term of order /3 2 we get, from the poles at z = — 2, 

r o Q ooo| (a)l = ( ' 7 ^4 (3)/32 ) Tr c (T a T b T c T d ) (fcf + 2k\ + fcf + 2k,k 2 + 2k 2 k 3 ) . (5) 

When the additional five permutations are added and the resulting four-point function inserted into Eq. (||), the 
contribution to the effective action from the quark loop becomes, 



r*| (a) (Ao) = 19 I drf (fx^A-A^Alf + {d^Mf - 2f a ™f bd ™(d i A a )-(d i A b )A c A%] , (6) 



after a partial integration. 

In the QED case we have N f Tr c (T a T b T c T d ) -> lj_and then our result for the two-derivative part of the QED 
effective action agrees with the earlier calculation in |19] . 



B. The pure Yang-Mills contribution 

We now turn to the pure Yang-Mills contribution, i.e. the diagrams (b)-(e) in Fig. 1. To evaluate the finite T part 
of these diagrams we will use the gauge condition d ■ A a = 0, and work in Feynman gauge. As in the fermion case, 
the functions gi = gi{y\ 1 ■ ■ . , fci, . . . , k 3 ) below are all functions of the relevant Feynman parameters and the external 
momenta. Their explicit forms can also be found in the appendix. Proceeding in a way similar to the previous section, 
we have for the ghost loop depicted in Fig. 1(b), 

r SoSo| (b)l = Tmiffff) / dyi / dy 2 / dy 3 / — r(z)C(z)cos[7rz/2](/3 2 ffl )^ /2 r((5-z)/2)r(z/2), (7) 

w 87r ' JO JO JO Jl+e-ioa z7r * 

where the color structure is //// = ff a j eb sfg ch f hd f = 5 ab 5 cd + S ad S bc + N(d abm d cdm - d acm d bdm + d adm d bcm )/4:, with 
N = 3 and dijk the completely symmetric structure constant. 
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For the graph in (c) we find, 



-^abcd\ I 0l 9 \ I I j„. / j. / j. / aZ 



«| (c)l = ) (////) J dm J q dm J q d V3 I -r mz) cob^- x 



<h 



- z/2 r(z/2) {5r((l-z)/2) + 16r((3-z)/2) + 32r((5-2)/2)} - g- {2+z)/2 T((2 + z)/2)x 



. 94 r((i-z)/2) - ^ 53 r((3-z)/2)l - | 953 - (4+2)/2 r((i-z)/2)r((4+z)/2) 



(8) 



where the color structure is as for the ghost loop. 
For the triangle diagram (d) we get, 



j^abcd I 
1 0000 



(d); 



= f (ff)(ff+ff) J o dyi J o dy 2 ^ ^- r ^ z ) cos[nz/2]((3 2 g 6 )-^ 2 x 

[9e\96 ~ 3. 97 )r((l~z)/2)r((2 + z)/2) - 15r((3-z)/2)r(z/2)] , (9) 



Where (//)(// + //) = fhfc fg df^a g efbhe + f ahmfb g mj = - 2 5 ab 5 cd - 5 ac 6 bd -8 ad 5 bc - Ndabmdcdm/2. 

Finally, we find for graph (e), 

/ 3?'n 4 \ f l i-l+e+ioo j 

rggoo| (e)l = (tAh (//+//)(//+//) / dl/i / — r(zK(z)cos[7rz/2](/3 2 .g 8 )- 2 / 2 r((l-z)/2)r(z/2) , (10) 

With (// + //)(// + //) = (p9e f bhe + f ahe f b 9 e ){f c g f f dhf + f chf f d gf) ^ 46 ab 5 cd + 26 ac 5 bd + 2S ad 6 bc + Nd abm d cdm . 

Contrary to the fermion case, at order O(/3 ) each diagram contains a logarithmic dependence on T and the external 
momenta kj. In addition, the T = part depends logarithmically on fcj and an ultraviolet cut-off A, that has to be 
introduced to regularize the loop-momentum integral. The logarithmic dependence on kj cancels out between the 
T = and T > parts for each permutation of each individual diagram, whereas the terms containing logT (log A) 
only cancel out in the total T > (T = 0) result, i.e. when all the different diagrams are added together. This was 
basically noted already in [{L8|| , and we have checked that it holds true in our calculations as well. There is also a 
linear divergence when (3k — » 0, at order (9(l//3), in all of the Eqs. (|7|)-(|l0|). This divergence originates from static 
propagators running in the loop, i.e. the propagators with a vanishing Matsubara frequency, oj n = 0. If only the 
non-static modes are integrated out, the static terms should be subtracted and our remaining result is then finite in 
the limit f3k — ► 0, as it should jL8|. We emphasize that the u> n = modes do not influence the two-derivative term. 

Taking into account all the permutations of Eqs. (|7|)-(^0|) and adding the different contributions, we find for the 
0(/3 2 ) term, 



Iffl-(e) (4>) = ~ 2567r 4 / dr\dx 



o 



2A%Al{d l A\) 2 + (dMo^o)) 2 + Yf acm f bdm ( d i A o)-(9iA b )A c A d 



(11) 



where we have performed an integration by parts. This result for the pure Yang-Mills contribution disagrees slightly 
with the previous finding from the background field method in |2^] , in that the first two term in Eq. (|ll]) are a factor 
(2/11) smaller, and the last a factor (— 1/2)Q. 

C. The total contribution 

By combining the results in Eqs. (||) and JTl|), the complete contribution to the effective action becomes, 

riVo) = r B | (o) (A ) + r E | (b) _ (e) (A a ) = 9 -^f - ij jfV J d^A-A^Alf + mA a M) 2 } - 



L E 



U9 J 3 qX (^f + 1 )l dT I dHf acm f bdm (d^y(dAW AZ] ■ (12) 



It should be noted, however, that the discrepancy is of marginal practical importance when it comes to the qualitative 
discussion of the influence of these QCD-terms in the Polyakov Loop action. 
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In the dimensionally reduced theory of QCD |l8|,[l9| , Eq. (|T^) provides the leading derivative interactions between the 
Ao-fields. How large the calculated derivative term is, compared to both the ones already present in Eq. (|^) as well 
as the omitted higher dimension operators, depends on the scales of interest. For instance, at the soft scale where 
di ~ gT and A) ~ T, we have g 4 2 d 2 A\ ~ g 6 T 4 , a factor g 2 higher than the term g 4 A\. 

In general, this effective theory is interesting on len gth scales \x\ <C (3 in the high temperature regime, especially 
when combined with nonperturbative lattice methods fl24|| . It is for example possible to study the non-perturbative 
Debye mass p5[ , and the 3d effective theory is also useful for calculations of the pressure in the quark-gluon phase, 
both perturbatively j^] and nonperturbatively pq| . Although the derived correction in Eq. (|l2| ) presumably gives a 
minor effect only, it is somewhat interesting to note that the terms are rather sensitive to the number of quark flavors. 



In Eq. (12), the coefficient of the first term is proportional to (7Nt — 9), which goes from —1 to 4/3 between Nf = 
and Nf = 3. Similarly, the second term in Eq. ( |l2] ) increases by more than 250% in the same range of Nf. This should 
be contrasted with the constant part of the ^-contribution to Eq. (||), that depends on Nf as (9 — Nf) and therefore 
only changes by 33% when going from Nf = to Nf = 3. Due to this strong behavior of the number of quark flavors, 
it is not inconceivable that the derivative interactions will make a small but noticeable difference between e.g. the 
pure glue theory and the three-flavor case. 

III. DERIVATIVE TERMS IN THE WILSON LINE MODEL 

The QCD dimensionally reduced theory describes accurately static phenomena at very high T, but the approxima- 
tions break down around a few times T c p4j . In addition, one is by construction omitting all dynamical information. 

To understand the features around T c a Ginzburg-Landau type of effective theory was proposed in [fl2|| . In this 
model, the potential is written in terms of I, 

V(l) = ai T 4 [-a 2 \l\ 2 -a 3 (l 3 + c.c.) + \l\ 4 ] . (13) 

The constants et^ are then used to fit the pressure above T c , with a 2 a function of temperature so that the global 
minimum of the potential is at I ^ (I = 0) above (below) T c [jl6| . One of the important aspects of the potential in Eq. 
( |l3| ) is the extremely rapid change around T c , due to a very sensitive dependence of a 2 on T/T c pq |. In a dynamical 
scenario one can therefore assume an instantaneous quench, where the value of I suddenly no longer corresponds to 
the correct minimum. The ^-field then rolls down the potential, and by coupling the Z-field to a linear sigma model 
the potential energy is converted into pions [ |l6| , |l7j ]. Even though the model is of phenomenological origin, it thus 
makes predictions that can be compared to experimental results. 

After the quench, the evolution of the Euler-Lagrange equations from the initial conditions requires, apart from 



the potential and the coupling to the chiral field, also a kinetic term for I [[16 17 1. Although the time dependence is 
beyond the calculation presented in this paper, we can provide the first perturbative QCD-correction to the spatial 
derivatives. The leading coefficient is the classical contribution to the derivative term, and comes from the kinetic 
term of A (x) in Eq. (||), as can be seen from the following argument ^7|: decomposing the Wilson Line in Eq. (|l|) 
into an octet L and the singlet I, 

L = L + 1— TrL = L+li , (14) 

where L is traceless and N = 3, we have 

Tr|9 4 L| 2 = Tr|9,L| 2 + 3|9^| 2 . (15) 

On the other hand, by a direct calculation in the static limit, 

Trl^LI 2 = g 2 {3 2 (diA%) 2 + Tr {commutator terms} . (16) 

When the commutator terms in Eq. (|l^) are rewritten in terms of L they can only involve the adjoint field, or products 
of I and L, since I (times the identity matrix) by itself commutes with all SU(3) matrices. Thus, by combining Eqs. 
( |l5| ) and (|l^), we have, 

(l/2)(d^ ) 2 - ^ [Tria.LI 2 + 3\dJ\ 2 ] + /(L) , (17) 

where /(L) corresponds to the commutator terms, rewritten as a function of L. At T c , g ~ 2.5 so that the leading 
coefficient for the kinetic term of / is 3/2g 2 ~ 0.3, which is reasonably close to the canonical value 1/2. Given the 
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unknown function /(L) it is not clear whether the kinetic term for the adjoint field actually is unique. However, L 
does not play any important role around T c , and can therefore be neglected on physical grounds fll^] . 

The procedure to obtain the kinetic term for I is thus to match terms in the effective theory T(^4o) to a corresponding 
\dil\ 2 piece in T(l). This means that the classical coefficient for \dil\ 2 will change when radiative corrections are taken 
into account in T(Aq). To lowest order, the kinetic term for Aq can receive corrections from the polarization tensor 
IP^P^I ) that would affect the <9J| 2 term via Eq. (17). However, with the optimal choice for the counterterms JT^] 
there is in fact no 0(g 2 )-contribution, so the renormalized kinetic term in T(Aq) remains (l/2)(diAo) 2 . 

Even though the kinetic term for Aq is unchanged at one-loop, this does not mean that \dil\ 2 is so. Since I = 
(f /3)Tr L, the |cV| 2 -term contains at least four powers of A$ when L is expanded in powers of A . The two-derivative 
term in T(l) can therefore receive a perturbative correction from the four-point function in Eq. (|^). Indeed, by using 
the relations 



AlA^diAlf 



--^| 2 -1] (Tr|a t L| 2 +3|^| 2 

144 W 



0{At) 



9 4 P 4 

2Trp,L, L]| 2 



0{A%) 
0{At) 



(18) 



we can rewrite Eq. (|2|) as 



(19) 



From the calculation of the four-point function we thus find the leading perturbative correction to the kinetic term 
of I. Including the contribution from Eq. (119), the kinetic term now becomes 



o T 2 



3T 2 
2g 2 



1 



g 2 m 

I6tt 4 



(7N, - 9) 



\dil\ 



3T 2 
2g 2 



\Sil\ 



(20) 



Since I is dimensionless, the correct dimension of the operator has to be supplied by some other scales. In the 
perturbative calculation the only scale is T, and hence the one-loop induced coefficient has the same T-dependence as 
the leading term in Eq. ([l?]). However, the perturbative correction does not depend on the QCD coupling constant 
and is therefore just a fixed number at this order. For example, for three flavors the coefficient is ~ 1.4 x 10~ 2 T 2 . 
Compared to the classical contribution (3/2g 2 )T 2 , the fraction of the one-loop correction is only O.Olg 2 . Even at T c , 
this is merely of the order of 5%, and at higher T even less due to the logarithmic decrease of g. 

Having derived the first correction to the kinetic term for I from perturbation theory, let us now discuss to what 
extent, and in what temperature range, the terms in Eq. (|l9| ) can be trusted. First of all, it should be noted that 
the form of the effective action in Eq. ([l9]) is not completely unique. The reason is that by a partial integration, 
and discarding any surface terms, we can always trade factors of [AgAg^Ap) 2 ] and (O^AqAq)) 2 for a term like 
[AgAgAg(<9 2 Ag)], but this equality does not hold at the level of I. In fact, if a term [AgAgAg(<9 2 Ag)] is kept in the 
action, not only do the coefficients in Eq. ( [l9| ) change, but there are also additional nonequivalent terms of the form 
(<9i|Z| 2 ) 2 and \l\ 2 [df(l + l*)]. Nevertheless, the action in Eq. ( |l9| ) is of course unique to order 0(A 4 ). Since higher order 
operators in the dimensionally reduced theory are further suppressed at high T, the predictions in Eq. ( |l9| ) should at 
the very least be reliable down to T > 2T C — 3T C , i.e. when the effective 3d theory itself is applicable. 

When T — > T c , the question is admittedly more subtle, as higher loop effects, higher dimensional operators and 
possibly nonperturbative effects become important. However, Eq. (|H)|) does not have to break down completely when 
the 3d theory does so. The 3d theory becomes invalid because the procedure of integrating out the non-static modes 
is unreliable when g 2 {T) T ~ ttT 0]. In contrast, T(l) is by construction valid near T c , so the question is rather how 
much the coefficient for the spatial derivative term changes. 

Considering first the operators of higher dimensionality, it is certainly possible to imagine that their bulk part 
follows from an expansion of Eq. (|l9|). The additional contributions that do not originate from these sources, e.g. 

terms like T 2 (|Z| 2 — l)"|<9iZ| 2 , that are at least of order A^ n+2 ^ (with n > 2), would then be suppressed. Not because 
they are unimportant a priori, but because their numerical coefficients are small. There are also higher derivative 
terms not accounted for in Eq. (19), like (|^| 2 — 1)(9 2 |/| 2 ) 2 , but they do not affect the kinetic term. To O(Aq), they 
are in fact straightforward to obtain from Eq. (0) and Eqs. (|7|)-(|To|). 
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When it comes to higher loop and nonperturbative effects, they will naturally induce a T-dependence in the radiative 
corrections to the kinetic term. Thus, (3T 2 /2 ff 2 )(l + c corr )|<9J| 2 -> (3T 2 /2g 2 )(l + c COTI (T/T c , T/A))\dJ\ 2 , which follows 
partly from the running of the QCD coupling constant in the higher loop effects. How much this will affect the kinetic 
term is difficult to estimate, but the small correction from the four-point function may indicate that the perturbative 
QCD-contributions are not too important for constructing T(l). 

IV. SUMMARY AND CONCLUSIONS 

In this paper we calculated the leading momentum dependence of the four-point function in QCD with Nf massless 
flavors, and related this contribution to terms in both the effective action T(Ao) and T(l). 

As for the derivative terms in r(^4o)j they will only have a minor influence when T 3> T c . As the temperature 
decreases and approaches T c the 3d theory becomes less reliable, but there could very well be a temperature region 
where the effective theory is still valid and the derivative interactions nonnegligible. Since the contribution has a 
rather strong dependence on Nf, a difference between the pure glue theory and e.g. Nf — 3 QCD could perhaps be 
noticed. 

From the four-point function, we also found the lowest one-loop QCD correction to the spatial derivative term in 
the effective theory T(l). The coefficient is independent of g and much smaller than the classical term, the ratio 
between the two being of the order of 10 -2 at T c . This derivation assumes that the strange quark mass m s can be 
neglected even at T c , which of course is an oversimplification, given that m s ~ T c . Nevertheless, the influence of m s 
is not likely to change the fact that the correction is small even at T c . 

At one-loop, there is an infinite number of terms that contribute to the coefficient of |9J| 2 , to the same order 
in g and T, as the four-point function. This follows from the fact the induced two-derivative interactions, with n 
external fields Aq 7 is of the functional form g n d 2 Aq /T^ l ~ 2 \ which corresponds to an expansion of I to at most order 
(n — 2) in the term T 2 \dil\ 2 . Some of these higher-dimensional contributions, maybe even the major parts, are already 
accounted for by rewriting the four-point function in T(Aq) in terms of I, as in Eq. (|l9|). In any case, since the 
four-point contribution is very small, it is reasonable to assume that the higher n-point functions give even smaller 
corrections. In that case, Eq. @ should give the correct order of magnitude for the total one-loop correction. 

As mentioned earlier, there are also higher loop effects that contribute to the kinetic term in T(l). For example, 
taking into account the two-loop correction to the diagrams in Fig. 1 gives c corr — ► c corr [l + ag 2 ]. To understand the 
reliability of the canonical term it is then crucial to know the magnitude of a. Surprisingly, studies of higher loop 
effects in the 3d effective theory r(^4 ) indicate that they only give corrections of the order of 30% at T c p4| . If 
these conjectures can be taken over to the Wilson Line model, one could in fact expect the derivative term (1/2)|9^| 2 
to change by perhaps at most a factor two, with all QCD-corrections taken into account. Of course, this has to be 
regarded as a highly speculative suggestion at the present stage. 

To complete the dynamical scenario one also needs the time dependence of I. Unfortunately, it is yet unclear how 
I generalizes to a real time formulation fl^ . Assuming that the form of the spatial derivatives can be extended to 
a Lorentz invariant form, the predictions for pion production and the evolution of |Z| will remain almost unchanged 
|^6|,0. In particular, if the Lorentz invariant kinetic term does not change by more than a factor of two, it can easily 
be compensated by a difference in e.g. the expansion rate of the plasma. 

Finally, to obtain a decisive estimate of the QCD-effects in T(l), one has to establish either a unique mapping from 
T(Aq) to r(i), or find a way to determine T(l) directly, perhaps numerically. Hopefully, the calculations presented in 
this paper can serve as a first step in that direction. 
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APPENDIX A: EVALUATION OF THE FEYNMAN DIAGRAMS 

In this appendix we outline our method for calculating the Fcynman diagrams shown in Fig. 1. We first follow the 
Feynman rules given in pq ] for Minkowski space-time. After all contractions and traces over spinor indices have been 
performed, we continue to a Euclidean space compact in the imaginary time direction: 
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Po 



where cj„ is the Matsubara frequency, u> n = (2n + 1)-kT (u> n — 2ttT) for fermions (bosons). 

Next, we use a Feynman parametrization to combine the denominators in the loop integral, and extract the T > 
part from the following relations p2|: 



T>0 



f-™t n B[f(Po) + f(-Po)} (bosons) 
/ST ^^[/(P0)+/(-P0)] (fermions) 



where ub = [exp(/3po) — 1] 1 , np = [exp(/3po) + 1] 1 an d 8 — ¥ + . The ghosts follow Bose statistics, despite their 
anticommuting properties. 

By shifting the vector momentum in the loop, p, to q = p + h(yi, ki), with h(yi, ki) a linear function of the external 
momenta ki, we then integrate over q. Finally, by performing a Mellin transform, 

rc±+ioo j 

{e*±l)~ l = —T(z)t(z)v ± x-' , (A3) 

Jc±—ioo ^ 7rl 

where v + = (1 — 2 1 ~ z ), u_ = 1, and with the contour specified by c + = e, c_ = 1 + e (where e — > + ), we can integrate 
over po after a change of variables. 

To illustrate the above procedure, consider the diagram in Fig. 1(e). Omitting the color factors for simplicity, and 
using the notation \p\ = p, \ki + hz\ — \k\2\ = *i2j we have to calculate the following integral: 

3<? 4 f d^p 1 1 T >o„. 4 f S+ZOC dp Q f 1 f d 3 p 1 

3ig 7T~ n B / "2/1 



(27r) 4 p2 _p2 p 2 _ {p+ J s _ too 2m J q J (27r )3 ^ 2 + {2 .. ^ + fe 2 2 } (1 _ yi) _ p 2]2 



1 pS+ioo /-oo 2 



3^4 ^ ,0 +tOO , . . , 

d yi / ^— n B / a^r^- 212 1 ( A4 ) 



2tt 2 7 y 5 _ ioo 2 ™ Jo [q 2 + 9%-pl? 

where we made t he s hift q = p + (1 — 2/i)*i2 and defined gs = yi(l — 2/i)*i2 in the last integral. After performing the 
q- integral in Eq. (A4), we are left with, 

3«5 4 f 1 , f s+loc d Po { 1 \ 1 3ig 4 f 1 , /■ 1+e+io ° , . 



8tt 7 7,5_ loo 2ttz Ve^-ly -p 2 8?r 2 J Q Ji +e - ioo 2?ri 7 ^/ 5s + 



3ig 4 f 1 J <fe 



167r 5 /2 7 n ./l +£ _ io o 2 ™ 



dl/i / — r(z)C(z)/3" 2 cos[^/2]. 9 - z/2 r((l - z)/2)T(z/2) , (A5) 



which is Eq. (10), except for the color factors. 

To check our method we also calculated the 0(Ji 2 k\ 2 ) contribution to the above diagram in a different way: we first 
performed the po-integral, by picking up the poles in the complex po-pkine, and then did the p-integral, without any 
Feynman parametrization, by expanding the integrand differently in different integration regions. The two results of 
course agree with each other. 

For completeness, we also give the functions /1, . . . , /j in Eq. (|j) and gi, . . . ,gs in Eqs. j7|)-(|To|). Using the conser- 
vation of momentum, and the shorthand notation ki + kj + ki — *«{, their explicit forms are as follows: 

fx = 2/2* 2 + 2/3*12 + (1 - y x - 2/2 - 2/3)^23 ~ [2/2*1 + 2/3*12 + (1 - 2/i - V2 ~ 2/3)^i23] 2 (A6) 
h = 3fc 2 + 4fcifc 2 + k\ + 2kik 3 + k 2 k 3 - 3(3*i + 2fc 2 + fc 3 )[y 2 fci + ^3^12 + (1 - Ift - 1ft - 2/3)^123] + 

+6[y 2 ki + y 3 ki2 + (1 - Vi - 2/2 - y3)ki2 3 } 2 (A7) 
h - (1/3) {5k 2 iyi (2 yi + k 2 2 [3 + 10yl - Wy 2 + 10y 2 2 + I0yi(2y 2 - 1)] + +5fc 3 2 [l + 2y\ + 2y 2 2 - 3y 3 + 2y 2 3 + 

+J/ 2 (4y 3 - 3) + yi(4y 2 + 4y 3 - 3)] + fc 2 fc 3 [8 + 20y 2 + 20y 2 - 10y 3 + 5y 2 (4y 3 - 5) + 5 yi {8y 2 + 4y 3 - 5)] + 

h[k 2 (3 + 20y 2 - 5y 2 + 5 yi (4y 2 - 3)) + fc 3 (2 + 20y 2 - 5y 2 - 5y 3 + 20 yi (y 2 + y 3 - 1))]} (A8) 
U = {[(1 - V2)k\ - 2/3*12 - (1 - yi - 2/2 - y 3 )ki 23 ] ■ [-1/2*1 + (1 - z)k 12 - (1 - yx -2/2 - y 3 )k 123 }} x 

{[-2/2*1 - 2/3^12 - (1 - 2/1 - 2/2 - 2/3)*m] • [~2/2*i - 2/3^12 + (2/1 + 2/2 + 2/3)^123]} - 

- {[-2/2*1 - 2/3&12 ~ (1 ~ 2/1 - 2/2 - 2/3)*123] ' [-2/2*1 + (1 - z)*12 ~ (1 ~ 2/1 ~ 2/2 ~ 2/3^123]} X 



8 



2/3)^123]} + 

- 2/2 - 2/3)^123]} x 

t- 2/3)^123]} 



(1 - 2/1 - 2/2 - 2/3)fc3] 2 - 



{[(1 - 2/2)^1 - 2/3^12 - (1 - 2/1 - 2/2 - 2/3)^123] • [-2/2&1 - 2/3^12 + (2/1 + 2/2 4 
+ {[-2/2^1 - 2/3^12 - (1 - 2/1 - 2/2 - 2/3)^123] • [(1 - Vi)ki - j/ 3 fci2 - (1 - 2/1 
{[-2/2% + (1 - z)fci2 - (1 - 2/1 - 2/2 - 2/3)^123] • [-2/2^1 - 2/3^12 + (2/1 + 2/2 - 
.91 = /l 

.92 = h 

g 3 = 3fc 2 + 2fe| + 2fc| + 3hk 2 + 3/ci/c 3 + 2k 2 k 3 + 2[(1 - y x )fei + (1 - 2/1 - 2/2)^2 
-(3fci + 2/c 2 + fc 3 )[(l - yi)h + (1 - 2/1 - 2/ 2 )fe 2 + (1 - 2/1 - 2/2 - 2/3)^3] 

.94 - (1/3) {kl(10yi - 20y 2 - 7) - 2fc 2 [10y 2 - 10y 2 + 10y 2 + 10yi(2y 2 - 1) - 3]- 
-fc 2 [17 + 20z/ 2 + 20y 2 - 30y 3 + 20jy 2 + 10y 2 (4y 3 - 3) + 10 yi (4y 2 + 4y 3 - 3)] - 
-2fc 2 fc 3 [2 + 20y 2 + 20y 2 - 10y 3 + 5y 2 (4y 3 - 5) + 5 yi (8y 2 + 4y 3 - 5)] - 
-2fc 1 [fc 2 (20j/ 2 - 52/2 + 5 yi (4y 2 - 3) - 3) + 5/c 3 (3 + 4y 2 - y 2 - y 3 + 4y 1 (y 2 + y 3 - 1))]} 

.95 = {[yifa + (1 + 2/1 + 2/2)^2 - (1 - 2/1 - 2/2 - 2/ 3 )fc 3 ] • [(2 - yi)h + (1 - y x - y 2 )k 2 + (1 - y x - y 2 - y 3 )k 3 }} x 
{\~y\k\ - (2/1 + j/ 2 )fc2 + (2 - 2/1 - 2/2 - 2/ 3 )fc 3 ] • [(2 - 2/1)^1 + (2 - 2/1 - 2/2)^2 + (2 - y\ - y 2 - y 3 )k 3 }} + 
{[(1 + yi)fci - (1 - 2/1 - 2/ 2 )fc 2 - (1 - 2/1 - 2/2 - 2/ 3 )fc 3 ] • [(1 + 2/i)fci + (1 + 2/1 + 2/ 2 )A: 2 + (1 + 2/1 + 2/2 + 2/3)^3]} x 
{[-2/ifci + (2 - 2/1 ~ 2/2)^2 + (1 - 2/i - 2/2 - 2/3)^3] ' [2/1^1 + (2/1 + 2/2)^2 + (1 + 2/1 + 2/2 + 2/3)^3]} (A14) 



(A9) 
(A10) 
(All) 

(A12) 



(A13) 



.96 = 2/1&12 + 2/2 fc 123 + (1 - 2/1 - 2/2)^ 2 234 ~ + fc 2 + (1 - 2/l)&3 + (1 - 2/1 - 2/2)fel] 5 

.97 = [(1 + 2/i)&3 - (1 - 2/i - 2/2)fc4][(2 - 2/i - y 2 )k 4 - yifo] , 



(A15) 
(A16) 



wher e, in order to simplify the permutations of the triangle graph, we did not use fc 4 = -{ki + fc 2 + fc 3 ) in Eqs. flAjg ) 
and (JA16D- Finally, 



,9s = 2/i (1 - 2/i)^i2 



(A17) 



See the homepage for the conference QM2001: http:/ /www, rhic. bnl. gov/ qm200\ 
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